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Docket No.: 740819-1036 



IN THE SPECIFICATION: 

Please amend the specification as follows: 

Please replace paragraph 1 on page 2 of the specification as follows: 

Specifically, the present invention provides a numerical calculation method, 
a numerical calculator or a recording medium that stores a numerical calculation 
program for a physical quantity U by solving A-U = f, wherein A is a coefficient 
matrix (in N rows by N columns; wherein N is a positive integer) obtained 
through discrete of a partial differential equation to be satisfied by the physical 
quantity U, and f is an inhomogeneous term (source term). In the numerical 
calculation method, the numerical calculator or the recording medium, the 
processes of setting 0 as an initial value of a number m of repeating times, giving 
0 as an initial value of a perturbation quantity (j> and setting (f - A.U°) as an 
initial value r° of a residual r; and repeatedly executing a first step and a second 
step while incrementing the number m of repeating times until an approximate 
solution U m is converged, and the first step includes the steps of setting an initial 
value U° of the physical quantity U; obtaining a predicted approximate value \|/ m 
of A- <j) = r m through repeated calculation performed by a first calculation unit 
including an internal solver, and the second step includes the steps of: obtaining, 
from the predicted approximate value \(/ m , a corrected approximate value <j> m for 
minimizing L 2 norm of a residual r^V^ 1 through an optimization routine 
performed by a second calculation unit; and giving (U m + <f> m ) as an approximate 
solution U m+1 and giving (r m - A* <f> m ) as a residual r m+1 , and in the second step, 
obtained elements of a vector sequence A- <f> m are sampled by a given sampling 
method to be stored in a memory, and a residual minimization coefficient a™ 
(wherein 1 = 1, . . ., L) used for obtaining the corrected approximate value <f> m is 
approximately obtained by using elements of a vector sequence A ^ k (wherein k 
= m - L + 1, . .., m - 1) stored in the memory. 
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Please replace paragraph 2 on page 5 of the specification as follows: 



Assuming that the number of repeating times is m, a synthesize 
perturbation quantity <p m for minimizing the L 2 norm of the residual i-™"*" 1 at 
(m+l)th-time and a new approximate solution U m+1 are defined as follows: 

<f> m =a"i// m +Y j a l m <f> m - M ... (5) 

wherein ai m (1 = 1, 2, 3, . . L) is a residual minimization coefficient and is a 
constant obtained through calculation described later. The residual is minimized 
so as to reduce the L 2 norm of the residual as r" 1 * 1 follows: 



Please replace last paragraph on page 8 through paragraph 1 on page 9 as 
follows: 



When this equation is expressed without using the symbol £, 

r m+l = r m _ (a m A gm + + a ™ A §™-2 + _ + a L m A8 m L+1 ) 

The residual minimization coefficient a™ can be obtained by using, for 
example, the least square method employing a condition for minimizing the 
(m+l)th-erdeF time residual norm [jr™^ 1 1| (a square root of a sum of squares). 
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When this equation is differentiated with respect to a™, the following is obtained: 
a/5cr / m (|r m+, f) = 0 
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